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We investigate the clumping of cold dark matter (CDM) at small scales. If the CDM particle is the neutralino, 
we find that collisional damping during its kinetic decoupling from the radiation fluid and free streaming introduce 
a small-scale cut-off in the primordial power spectrum of CDM. This cut-off sets the scale for the very first CDM 
objects in the Universe, which we expect to have a mass of ~ 1O" 12 M . For non-thermal CDM candidates, such 
as axions, wimpzillas, or primordial black holes, the cosmological QCD transition might induce features in the 
primordial spectrum at similar mass scales. 



1. Introduction 

Dark matter with an equation of state that has 
been non-relativistic (p <C p) since structure for- 
mation started (around matter-radiation equal- 
ity) is called cold dark matter 0. CDM gives rise 
to the hierarchical formation of large structures 
(> 1 Mpc) in the Universe, i.e. small structures 
form first and grow to larger structures later. 
The growth of CDM density fluctuations is sup- 
pressed during the radiation dominated epoch. 
Thus the mis CDM density fluctuations go like 
k 3 at large scales, but increase only logarithmi- 
cally with wavenumber at scales well below the 
Hubble scale at matter-radiation equality. 

At very small scales (<C 1 Mpc), the power 
spectrum and the evolution of CDM density fluc- 
tuations has not been discussed in detail so far, 
although the understanding of the small-scale be- 
havior of CDM density fluctuations is essential 
for a realistic estimate of expected rates in CDM 
searches. The reason is that analytic calculations 
in this deeply non- linear regime are very challeng- 
ing and numerical simulations do not have the 
dynamical range to resolve scales as small as the 
solar system. A priori one can say that there has 
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to be a cut-off in the CDM power spectrum at 
some very small scale, otherwise the energy den- 
sity in the fluctuations itself would be infinity. 

The nature of CDM is unknown. Here we men- 
tion three popular candidates with very differ- 
ent properties: The first one is the lightest neu- 
tralino, which probably is the lightest supersym- 
metric particle j2|. Its mass is expected to be in 
the range 40 GeV - 600 GeV §. The neutralino 
is a mixture of the neutral gauginos and higgsinos, 
thus it interacts through weak interactions only. 
Below we assume that the lightest neutralino is 
the bino, because in the constrained minimal su- 
persymmetric standard model the dominant con- 
tribution in the mix comes from the bino J3j . 

A second CDM candidate is the axion ff The 
axion mass has been restricted to the be 10 -6 eV 
- 10~ 2 eV HQ and contributes to the dark matter 
if the mass is small. Axions interact much weaker 
than weakly interacting particles, the interactions 
are suppressed by the Peccei-Quinn scale, which 
is ~ 10 12 GeV. Thus axions are never in ther- 
mal equilibrium with the radiation fluid. An ex- 
ample of small-scale structure in CDM has been 
found from the initial misalignment mechanism 
of axions by Hogan and Rees ||. It turns out 
that, if the Peccei-Quinn scale is below the reheat- 
ing temperature after inflation, large isocurvature 
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perturbations in the axion density are created 
once the axion mass is switched on during the 
QCD transition. It has been shown that ax- 
ion mini-clusters with ~ 10 _12 Mq and radii of 
~ 0.1-R© might emerge, which may be observed 
by means of pico- and femtolensing 0. 

Let us mention a third CDM candidate: pri- 
mordial black holes Their mass should be 
> 10~ 16 Af Q in order to survive until today ||]. 
Primordial black holes interact with the rest of 
the Universe via gravity only. They may be found 
or excluded by gravitational lensing Jl^Q . 

2. Damping scales for CDM 

For any thermal CDM species there are two 
mechanisms that contribute to the damping of 
small-scale fluctuations: During the kinetic de- 
coupling of CDM from the radiation fluid the 
mean free path is finite and thus collisional damp- 
ing occurs. After the kinetic decoupling has been 
completed free streaming can further wash out 
the remaining fluctuations. 

For neutralinos, chemical freeze out happens 
at ~ m/20 > 2 GeV In contrast kinetic de- 
coupling happens at much smaller temperatures 
fPH , because elastic interactions with the radia- 
tion fluid are possible at temperatures as low as 
1 MeV. However, due to the large momentum of 
the neutralinos, many collisions are needed for a 
significant change of the momentum. It turns out 
that N ~ m/T collisions can keep the neutrali- 
nos in kinetic equilibrium and thus the relaxation 
time can be estimated as r ~ iW co u, where the 
collision time for a bino is given by 
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with M being the slepton mass and m the mass of 
the bino. With a slepton mass M = 200 GeV and 
a bino mass of m = 100 GeV, the relaxation time 
is given by r ~ (10 MeV/T) 4 tn, where in is the 
Hubble time. Kinetic decoupling of neutralinos 
happens at T ~ 10 MeV. 

We incorporate dissipative phenomena by de- 
scribing the CDM as an imperfect fluid |l2| . The 
coefficients of heat conduction, shear and bulk 
viscosity are estimated to be mx ~ rj ~ £ ~ nTr, 



where n is the number density of CDM particles. 
We find that the damping of density perturba- 
tions goes as 
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where the damping scale depends on the mass 
of the neutralino and the slepton masses and is 
typically M D = 10" 13 M Q - 10- 10 M Q . For com- 
parison, the CDM mass within a Hubble patch is 
- 10- 4 M Q at T ~ 10 MeV. 

Free streaming leads to additional damping. 
The velocity of neutralinos right after kinetic de- 
coupling is v ~ (T/m) 1 / 2 ~ 10~ 2 . Free stream- 
ing also gives rise to exponential damping, due to 
the velocity dispersion. The typical free steam- 
ing scale is 10" 12 M Q - 10- 10 M o . Thus both 
damping mechanisms operate approximately at 
the same scale. The power spectrum of neutralino 
CDM is cut off at M < 10~ 12 M Q . 

The mechanism of collisional damping also 
works for CDM in the form of a heavy neutrino 
(~ 1 TeV), but it does not work for wimpzillas 
fi3| |, because these are too heavy to ever be in 
thermal equilibrium. Free streaming induces a 
cut-off in the power spectrum for all mentioned 
CDM candidates. The scale and the strength of 
the damping depend on the masses and the pri- 
mordial velocity distributions. 



3. QCD induced CDM clumps 

Besides damping mechanisms there are also 
processes that might enhance the primordial 
CDM spectrum at small scales. One is the forma- 
tion of axion mini-clusters jf| that we mentioned 
already in the introduction. 

Together with Schmid and Widerin one of the 
present authors has found that large amplifica- 
tions of density fluctuations might be induced by 
the QCD transition at scales 10~ 20 M© < M < 
10- 10 M Q @@- The mechanism is the fol- 
lowing: During a first-order QCD transition the 
sound speed vanishes. Thus the density pertur- 
bations of the dominant radiation fluid go into 
free fall and create large peaks and dips in the 
spectrum, which grow at most linearly with the 
wavenumber. These peaks and dips produce huge 
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gravitational potentials. CDM falls into these 
gravitational wells. It is important to note that 
this amplification mechanism works for matter 
that is kinetically decoupled at the QCD tran- 
sition around 150 MeV. For neutralinos this is 
not the case. Large inhomogeneities in the neu- 
tralinos would be washed out by collisional damp- 
ing later on. A structure similar to the acoustic 
peaks in the photon-baryon fluid might survive. 
The large inhomogeneities in the radiation fluid 
are completely washed out during neutrino de- 
coupling at ~ 1 MeV. 

4. The first CDM objects 

The smallest scales that survive damping are 
the first scales that go non-linear, thus these 
scales form the first gravitationally bound objects 
(apart from primordial black holes, if they exist) 
in the Universe. Let us estimate their size if the 
CDM is the neutralino. The scale is given by the 
cut-off at M ~ 1(T 12 M . With a COBE normal- 
ized |l5| CDM spectrum we find the rms density 
fluctuations at equality to be 
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which is 2 x 10 _2 (10~ 1 ) for the spectral index 
n = 1(1-2). Thus these objects go nonlinear at or 
shortly after equality, at a redshift of ~ 10 2 (10 3 ). 
If we assume that these clouds are spherical they 
would collapse to a radius of 10 3 (10 2 )i? Q today, 
which is by chance a very interesting scale for 
observations. Although 10 _12 A/ Q in a volume of 
(10 2 i? Q ) 3 seems to be an extremely diluted cloud, 
the overdensity of such a cloud today would be 
~ 10 11 . In more optimistic scenarios (larger tilt 
and/or additional peaks in the spectrum from the 
QCD transition) it is even possible that these 
clouds are so compact that pico- and femtolensing 
can be used to search for them. 

It is unclear whether some of these first objects 
can survive up to today, this will be the subject 
of further studies. We think that understanding 
and revealing the small-scale structure of CDM 
might help us to learn more about the nature of 
CDM. 
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